Theory of the Ordered Phase in A-site Antiferromagnetic Spinels 



SungBin Lee 

Department of Physics, University of California, Santa Barbara, CA 93106-9530 

Leon Balents 

Kavli Insitute for Theoretical Physics, University of California, Santa Barbara, CA 93106-9530 

(Dated: August 21, 2008) 

Insulating spinel materials, with the chemical formula AB2X4, behave as diamond lattice antiferromagnets 
when only the A-site atom is magnetic. Many exhibit classic signatures of frustration, induced not geometri- 
cally but by competing first and second neighbor exchange interactions. In this paper, we further develop a 
theory^ of the magnetism of these materials, focusing on the physics observable within the ordered state. We 
derive a phenomenological Landau theory that predicts the orientation of the spins within incommensurate spi- 
ral ordered states. It also describes how the spins reorient in a magnetic field, and how they may undergo a 
low temperature "lock-in" transition to a commensurate state. We discuss microscopic mechanisms for these 
magnetic anisotropy effects. The reduction of the ordered moment by quantum fluctuations is shown to be 
enhanced due to frustration. Our results are compared to experiments on MnSc2S4 , the best characterized of 
such A-site spinels, and more general implications are discussed. One prediction is that magnetically-induced 
ferroelectricity is generic in these materials, and a detailed description of the relation of the electric polarization 
to the magnetism is given. 



I. INTRODUCTION 

Frustrated magnets, in which competing exchange interac- 
tions cannot be simultaneously minimized, have long been a 
subject of theoretical and experimental study. Their funda- 
mental interest comes from their tendency to show more pro- 
nounced effects of fluctuations than their unfrustrated counter- 
parts, and from prospects of observing exotic ground states as 
a consequence of frustration-induced sensitivity to weak per- 
turbations. From a more practical materials science perspec- 
tive, they are of particular recent interest because they tend to 
display non-collinear magnetic ordering. Such non-collinear 
ordering is quite generally connected to magnetically-induced 
ferroelectricity, making frustrated magnets a rich and produc- 
tive hunting ground for multiferroics. In this paper, we study 
a particular class of frustrated spinel materials, with the chem- 
ical formula AB 2 X 4 , in which only the A atom is magnetic. 
Such materials are described as antiferromagnets on the dia- 
mond lattice. Somewhat surprisingly, although the diamond 
lattice is not geometrically frustrated and admits a simple two 
sublattice collinear Neel state, many of these A-site magnetic 
spinels do exhibit significant signs of frustration. This in- 
cludes a large ratio ("frustration parameter") / = \Qcw\/T c 2 
between the Curie-Weiss temperature <dcw an d an ordering 
or freezing temperature T c . For example, experiments find 
/ w 10 - 20 in C0AI2O4™, and / w 12 in MnSc 2 S 4 ! 5 
A recent theoretical study attributes this to the competition 
between first and second neighbor exchange interactions, J\ 
and J 2 , which can be comparable in these materials!^ The- 
oretically, for J2/J1 > 1/8, the classical ground state be- 
comes highly degenerate, consisting of coplanar spirals whose 
wavevector can be arbitrarily chosen on some "spiral surface" 
in momentum space. This degeneracy was suggested to be 
responsible for the observed signs of frustration, including 
large /, prominent diffuse neutron scattering in the param- 
agnetic state, and some low temperature specific heat anoma- 
lies. While encouraging, many of the predictions of this the- 



ory cannot currently be tested due to the absence of single 
crystal neutron scattering data. 

In this paper, we develop this theory further, in order both 
to capture more detailed physical properties of this class of 
materials, and to make further predictions which might more 
readily be compared to existing and future experiments. We 
focus on physics than can be directly observed in the ordered 
state, which has been fairly well characterized in MnSc2S4 . 
Specifically, we consider details of the magnetic anisotropy, 
and the magnitude of the local ordered moments at low tem- 
perature. The theory of Ref. Q] was based on an Heisenberg 
model, which possesses 0(3) (or SU(2)) spin rotation sym- 
metry and hence exhibits no preference for the absolute ori- 
entations of the spins themselves in the ordered state. Experi- 
mentally, in MnSc2S4 the spins are observed to lie in a definite 
plane. Moreover, the ordering wavevector describing the axis 
and pitch of the spiral in real space displays a "lock-in" behav- 
ior at low temperature, in which it becomes commensurate 
with the underlying spinel lattice. In the Heisenberg model, 
there is no explanation for this lock-in. We show here that both 
the choice of spiral plane and the commensurate lock-in of the 
spiral wavevector can be understood by considering magnetic 
anisotropy effects. By an extended phenomenological Landau 
analysis, we can describe the magnetic orientation selection 
across the broader family of A-site spinels - which has not 
yet been studied experimentally - and predict some interest- 
ing "spin flop" and reorientation effects in applied magnetic 
fields. We also consider, as mentioned, the value of the or- 
dered moment, which experimentally shows a relatively large 
(for an S = 5/2 spin) 20% suppression from the classical 
value. We show that, despite the large Mn 2+ spins, this can ac- 
tually be accounted for by quantum fluctuations, provided fur- 
ther neighbor interactions are sufficiently small, due to the en- 
hancement of fluctuations by frustration. Finally, we discuss 
the microscopic mechanisms behind the magnetic anisotropy 
of these materials, which may arise both from dipolar interac- 
tions and spin-orbit effects. In MnSc2S4 , we find that spin- 
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FIG. 1: The diamond lattice with the 1st, 2nd and 3rd nearest- 
neighbors coupling Ji , J2 , J3 respectively. 



orbit induced exchange anisotropy is the only one of these two 
mechanisms consistent with experimental observations. 

We emphasize that though we pay particular attention to 
the comparison with MnSc2S4 , the A-site spinels comprise a 
quite large set of interesting magnetic materials, and the the- 
oretical analysis of this paper is formulated in such a way as 
to apply to the entire family. It therefore has numerous im- 
plications for many materials, and should be quite useful as a 
guide to future experiments. Of particular interest is the possi- 
bility of observing ferroelectricity and magnetoelectric effects 
in these compounds. Our modeling of magnetic anisotropy 
contains the essential ingredients for a theory of magnetically- 
induced ferroelectricity. We present some basic observations 
of this type in the Discussion at the end of the paper. 

The remainder of the paper is organized as follows. In 
Sec.[nJ we describe a phenomenological form of the magnetic 
anisotropy in terms of the order parameter, based on symme- 
try constraints, and the resulting ground states. In Sec. Ill we 



discuss the magnetization process and a spin flop transition in 
a field. Sec. IV discusses tendency of the spiral wavevector 
to lock to commensurate values, and associated phase tran- 
sitions. We show in Sec.[V|how quantum fluctuations can be 
included in the theory. Then, in Sec.[Vl]we consider the possi- 
ble microscopic sources of the magnetic anisotropy, and con- 
clude that in M11SC2S4 , it is most likely dominated by spin- 



orbit induced exchange anisotropy. We conclude in Sec. VII 
with a summary of results, and a discussion of experimental 
phenomena, including magnetically induced ferroelectricity. 
Some technical calculations are included in the Appendices. 



II. SPIRAL SPIN STATE AND SPIN ROTATIONAL 
SYMMETRY BREAKING 

A. Heisenberg model and its ground states 

A minimal Heisenberg model description for the mag- 
netism of these materials was studied in Ref.[U Here the spins 
reside at the spinel A sites, which form a diamond lattice(see 
FigjTJi, composed of the 2 interpenetrating fee lattice. The 
Hamiltonian, in zero magnetic field, is simply 



-^hcis 2 / j Jij^i ' ^J* 



(1) 



Here we consider classical unit vector spins |Si| = 1. We 
consider coupling between up to third neighbor diamond sites, 
i.e. Jij = Ji, J2, J3 for first, second, and third neighbor sites, 
respectively. Though the diamond lattice with only nearest 
neighbor spin exchange J\ has an unfmstrated unique ground 
state, the inclusion of additional interactions (2nd,3rd near- 
est neighbor etc) rapidly produces frustration. Following the 
logic of Ref. 1 we presume that the first second nearest neigh- 
bor exchanges, Ji , J2, are dominant, and treat the third neigh- 
bor coupling J 3 as a small (but important) degeneracy break- 
ing perturbation. 

Ground states of this Hamiltonian can be found for arbitrary 
Ji by the method of Luttinger and Tisza. They take the form 1 
of coplanar spirals 



sf B) = Ide^^ 2 + c.c. 
2 



(2) 



where the order parameter d is a complex 3-component vector 
satisfying 



d-d = 0, 
d-d* = 2. 



(3) 



These two constraints, and the choice of 7, ensure that the 
magnitude of each spin is unity, |S,| = 1. One has 



7 = Arg 



(4) 



where the sum ^ 1S taken over sites is taken over all sites j 
on the B sublattice, with i fixed as an arbitrary A sublattice 
site. The physical meaning of d is made clear by solving the 
constraints: 



d = ei + ie 2 , 



and defining 



e 3 = ei x e 2 



-d x d* 



(5) 



(6) 



Here e 1; e2,e3 are three mutually orthogonal unit vectors. 
The first two span the plane in which the spins reside, and 
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§3 is the unique normal to the plane. A phase rotation of d ro- 
tates the spins within the plane, or equivalently translates the 
spiral along its axis, while leaving the spin plane and hence e 3 
unchanged. 

The energy of spiral states of this type is readily evaluated. 
It is sufficient to linearize in J3, in which case one finds the 
energy per unit cell (this is twice the energy per spin) 



Ej(k) = E 12 (k) + E 3 (k), 



(7) 



exchanges and the smaller J3 exchange, respectively. Explic- 
itly, 



E12 
SE 3 



= J: 



16J 2 (A(k) 
A(k)' 
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8J 2 



4J 2 
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(8) 
(9) 



where En and E3 are the contributions from the large J\ , J 2 
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(10) 



Treating J 3 perturbatively, we first minimize i? 12 . For 
Jz/\Ji\ < 1/8, the minimum occurs for k = 0, while for 
Jz/\Ji\ > 1/8, it occurs along the surface defined by A(k) = 
I «7i |/8^2- In the latter case, the 3rd nearest-neighbor exchange 
breaks the "spiral surface" degeneracy. A combination of ana- 
lytical and numerical arguments (see Appendix [A) determine 
the selected wavevectors on the spiral surface. We assume 
antiferromagnetic J3 > 0, in which case the minimum en- 
ergy is realized with a wave vector of the form q = (q, q, k), 
where the relation of k to q varies depending upon the magni- 
tude of J2/ J\. The direction of the wavevector thereby varies 
from the (111) to the (110) directions, with an intermediate 
(111*) region in which the k is chosen as close as possible to 
q, since the (111) directions do not intersect the spiral surface. 
See Fig |2] and Appendix |A| for further details. We note that 
this wavevector, determined from the third nearest-neighbor 
exchange J3, is different with the one determined by thermal 
fluctuations.^ 

For the specific material MnSc2S4 , the magnetic structure 
is known from neutron diffraction. 6 At low temperature the 
ordering wavevector is k = q = 37r/2(l, 1,0), and the re- 
finement indicates ferromagnetic Ji < 0. Comparison to the 
theoretical structure and the measured Curie- Weiss tempera- 
ture allows one to constraint the couplings 1 . When J3 is very 
small, one has J\ ~ — 10.5K and J2 ~ 8.75K. More gener- 
ally, fixing k = 37r/2(l, 1, 0), one has 



J*/\Jl\- — i ^ — ■ (11) 



From this relation, Ja/|Ji| varies from 0.88 to 0.94 when 
Ja/\Ji\ is increased from 0.01 to 0.04. 
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FIG. 2: The selected wavevector of the diamond antiferromagnet for 
antiferromagnetic J3. We plot q/iv as a function of J2/J1, where 
the ground state wavevector has the form (q,q,k). The direction 
(choice of k) is indicated by the labels (111), (111*), (110)- see text 
and Appendix |7£] for details - in each of the regions separated by 
vertical lines. The first four lowest order commensurate wavevectors 
q m ,n for which lock-in transitions are expected are also indicated by 
labeled dots (see Sec.jlVj. 



B. Magnetic anisotropy 

The Heisenberg model leaves the plane and phase of the 
spin spiral undetermined, because they can be continuously 
rotated using the SU(2) symmetry of the Hamiltonian. In re- 
ality, this symmetry is broken by the crystal lattice and "spin 
orbit" effects (in fact arising both from quantum mechanical 
spin-orbit coupling and dipolar interactions between spins) 
that couple spin and spatial rotations. Indeed, in MnSc2S4 , it 
is known that the spins in the (110) spiral lie in a (001) plane. 
This is determined by physics outside the Heisenberg model. 
Furthermore, the commensurate magnitude of the wavevector 



- q = go(l) 1) Q) with qo = 3ir/2 exactly within experimen- 
tal resolution - is also related to anisotropy effects. In the 
Heisenberg model, obtaining this value of q Q at T = re- 
quires fine-tuning of the ratio of J2/J1, and even with such 
tuning, the magnitude would generally deviate at T > 0. 

To understand these effects, we first adopt a phenomenolog- 
ical Landau theoretic approach constrained only by symme- 
try. This consists of time-reversal invariance, which reverses 
spins, and the space group, Fd3m, of the spinel lattice. The 
full space group is generated by 6 operations, which may be 
expressed in terms of translations T t by the vector t, rotations 
lZ n [8] by angle 8 about the n axis, and the inversion X about 
the origin. In our coordinate system, the generators Qi are 
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(14) 

(15) 
(16) 
(17) 



Because the spin is a pseudo-vector, its transformation under 
each of these operations is given by 



S(x) -> Det[6] ■ 6 1 ■ S(6 ■ r + 1) 



(18) 



where O is the orthogonal matrix giving the rotation/inversion 
part of the operation (r — > O r) and t is the translation vector. 

We are interested in the effect of spin-orbit coupling within 
the ordered phase of these materials. In this case, the sym- 
metry is already reduced from that of the full crystal by the 
magnetic order. Specifically, we assume an ordered state of 
the form predicted by the Heisenberg model, i.e. satisfying 
Eq. |2]) with k determined to be one of the values selected by 
Ji,Ji,Jz> but with d arbitrary up to the constraints in Eq. p). 
We seek a Landau free energy as a function of d. Since we 
restrict to a fixed k, we should consider only those symmetry 
operations which leave k invariant (up to inversion). This is 
the little group of the wavevector k. Under each element in 
this little group, because the wavevector is invariant, one can 
define a corresponding transformation for d, under which the 
free energy must be invariant. 

We consider the two major regimes of phase space in which 
the form of k is simple. For 1/8 < J2/I Ji| < 1/4, we have 
k = fc(l, 1, 1). The little group is generated by the transfor- 
mations 03,05, Qq in this case. Under these operations, the 
order parameter transforms according to 



03 : di 

05 :d 

06 :d 



d 3 , 
d*, 
e ik d. 



di, 



(19) 



In the case J2/J1 > 0.7, one has k = fc(l, 1,0), for which 
the little group is generated instead by Q\, Q4, 05, Qq. Under 
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FIG. 3: Directions of the normal §3 to the plane of spin ordering 
selected by magnetic anisotropy terms in the cases (a) of a (111) 
wavevector and (b) of a (1 10) wavevector. In (a), the symbol _L( 1 11) 
indicates that any plane with S3 • (111) = is a ground state. 



these operations, we find 
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(20) 



Using these symmetries, we can determine the most general 
allowed form of the free energy at any given order in d, for 
each of these two cases. Our focus is on terms which violate 
SU(2) symmetry, induced by spin-orbit coupling or dipolar 
interactions. As usual within Landau theory, we expect terms 
which involve smaller powers of the order parameter to be 
most important. We therefore consider the leading quadratic 
terms other than the trivial |d| 2 one. For the k = (k, k, k) 
states, we find a single non-trivial invariant: 

/ m (d) = c + d 2 ) + d* 2 (d ± + d 3 ) + di(d 2 + d 3 )} . 

(21) 

For the wave vector k( 1,1,0), the quadratic free energy con- 
tains two non-trivial invariants: 



/110(d) = ci(d\d 2 + c.c) + c 2 dld 3 . 



(22) 



These quadratic terms distinguish different planes in which 
the spins spiral energetically. We note that both /m and /no 
are invariant under arbitrary phase rotations of the d fields. 
Physically, this implies rotations of the vectors ei and 62 
within the plane normal to §3 cost no energy. Therefore we 
expect that these terms may be rewritten in terms of 63 alone. 
This is indeed the case. To do so, it is convenient to intro- 
duce a parametrization of d which solves the constraints in 
Eqs. Q: 



"737 ^7 



(23) 
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where we have defined the spinor z a , 

z = (e* cos0,e^ 2 ; 



(24) 



which satisfies \z\ | 2 + \ zi | 2 = 1. Here er is the vector of Pauli 
matrices, and e a p is the anti-symmetric matrix with ei2 = 1. 
It is straightforward to show that 



e 3 = ZaVafSZp. 



(25) 



By explicit evaluation using Eqs. ( 23|24 l, one can readily 
show 



/m = c l-[ef + e | + e 3 T 



/ 11 o = -2 Cl e^ + C2 [(e^) 2 + (e!) 2 l 



(26) 
(27) 



Now the energetically preferred plane for the spins is ap- 
parent. They are illustrated in Fig. [5] For k = [k, k, k), 
the ground state has 63 = (1, l,l)/v3 for c > 0, and 
§3 ■ (1, 1,1) = for c < (i.e. in the latter case, the vec- 
tor e 3 is still free to rotate anywhere within a plane). For 
k = (k, k, 0), three distinct directions of §3 are possible de- 
pending upon the values of c\ , c 2 - see Fig. |3]for details. 

At this stage it is possible to compare with experimental re- 
sults on MnSc2S4 . Refined neutron scattering data in Ref. |6 1 
indicated spiral order of the type discussed here with wavevec- 
tor q = (q, q, 0) and spins aligned within the (001) plane. We 
see that the Landau theory indeed captures this order, provided 
the phenomenological parameters c\ , c 2 are taken to lie within 
region I of the phase diagram in Fig. [3] Note that this is not 
"fine-tuning", as this region occupies a finite fraction of the 
phase diagram. However, it is still interesting to understand 
the microscopic reason for the system to be in region I rather 



than II or III. We will return to this question in Sec VI 



III. MAGNETIZATION PROCESS 

In this section, we consider the evolution of the spin state 
in an applied magnetic field. Neglecting magnetic anisotropy, 
we may expect a smooth evolution, in which the spins adopt 
a canted (conical) configuration with a non-vanishing compo- 
nent along the field, and this canting gradually increases until 
the spins become fully aligned at saturation. In the presence 
of magnetic anisotropy, however, the spins have an intrinsic 
preference for particular planes, which, in some field orienta- 
tions, competes with the tendency of the spins to adapt to the 
field. We study these two situations below. 



A. Heisenberg model 

We first neglect magnetic anisotropy and consider simply 
the classical Heisenberg Hamiltonian with an added Zeeman 
magnetic field 



H 



J,h 



(28) 



h s J(J\+J 2 ) 




hKh+h) 



0.1 0.2 0.3 0.4 0.5 



FIG. 4: Saturation fields h aa t/(Ji + J2) as a function of J2K J1 + 
J 2), for the ferromagnetic (lower curve) and antiferromagnetic (up- 
per curve) cases. 



We seek ground states with normalized spins | S; 
following ansatz: 



S A(B) = Ide i(k ' Xl±7/2) 



c.c + m, 



with the constraints 



d d 
d ■ m 



0. 
0. 



d • d* 



m 2 = 1. 



1, using 
(29) 

(30) 
(31) 

(32) 



We now evaluate the energy for these states. It is necessary to 
consider ferromagnetic and antiferromagnetic J\ separately. 



1. Ferromagnetic Ji 

In the ferromagnetic case, evaluating the energy per unit 
cell using the Hamiltonian in Eq. d28l, one obtains 



1 



E$£ = -E,,(k)\d\ 



\m\ 2 Ej(0) - 2h • m. 



(33) 



Here Ej(k) is the energy function for a pure spiral in zero 
field, given in Eq. dTJ. 

This energy function is minimized as follows. Only the 
third term is dependent upon the orientation of d and m, and it 
is minimized if we choose m = mh along the field direction. 
Then we must choose, similarly to Eq. Q, 



d = \J\ - m 2 (ei + ie 2 ), 



(34) 



with e 3 = {?! x e 2 = h. This indeed describes a conical spin 
state. For fixed m and hence |d| 2 = 2(1 — m 2 ), the energy is 
minimized by the wavevector k = q which minimized Ej(k). 
This implies that the wavevector is independent of magnetic 
field. Finally, we can minimize over m, which gives 



h 

h S at 



(35) 



which is valid for fields below the saturation field, which in 
this ferromagnetic case is 



h FM 



Ej(0) - Ej(q) = AE. 



(36) 
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Here we define AE for later convenience. We see that the 
magnetization increases perfectly linearly up to saturation. 
The saturation field itself varies with the exchange couplings 
and in particular J 2 / J\ in a non-trivial manner as the ordering 
wavevector q varies - see Fig. [4] Since the ground state itself 
is ferromagnetic for J 2 < J\) 8, the saturation field vanishes 
in this region. 



2. Antiferromagnetic Ji 

Next consider the case of antiferromagnetic J\ . In this case, 
the energy function is 



E 



AFM 



1 



Jih = -Ej(k)\d\ 



-(X/(0) + 8Ji)|m| 2 -2h-m. (37) 



The difference from Eq. ( |33| > can be understood as arising be- 
cause of the cost 8J1 of flipping the four nearest-neighbor 
bonds per site from anti-parallel to parallel spin alignment. 
Repeating the analysis of the previous subsubsection, we 
again find a linear magnetization curve (i.e. Eq. (35 1), but 
with the saturation field 



AFM 



8 Ji + AE. 



(38) 



B. Anisotropy and spin flop transition 

We now turn to the effects of magnetic anisotropy, and in 
particular the competition between the magnetic field and the 
intrinsic preference for the spin ordering plane. Lacking a 
microscopic model for the anisotropy, we cannot reliably ex- 
plore the full phase diagram for all fields. However, since we 
expect that the anisotropy is relatively weak compared to the 
exchange, the portion of phase space in which the field and 
anisotropy are actually competitive is restricted to small fields. 
In this regime, the contribution of the anisotropy to the energy 
should be approximately unchanged from that at zero field, 
and hence we may model it by the same phenomenological 
function given in Sec. |I1 B| That is, we add to the Heisenberg 
energy Ej^ the terms /m, /no, as appropriate. 

Let us focus on ferromagnetic J\ with q = (q, q, 0) for 
simplicity. The discussion is not significantly modified in the 
antiferromagnetic case. The energy function is now 



E, 



™ =^(k)|d| : 



|m| 2 i; J (0)-2h-m + / 11 o[d]. (39) 



In small fields, we may fix k = q the zero-field ordering 



wavevector which minimized Ej. We can use Eq. (34 1, with 



however §3 = rh not necessarily parallel to h. Inserting this 
into the energy, we find 



E, 



FAI 



= Ej(q) + AE\m\ 2 2hm 
+/no[e 3 = rh]. 



(40) 



where /no [63] is given in Eq. (27 1. Here we have approxi- 
mated m ss in the anisotropy term, since the neglected cor- 
rections are of 0(m 2 ci. 2 ), i-e. small both in the magnetization 
and the anisotropy. 



We can now minimize Eq. (40 » over the magnitude of the 
magnetization at fixed orientation, which gives 



(41) 



AE ' 



and the energy, which now depends only upon the orientation 
rh: 



(42) 



up to constants independent of rh. We caution that in these 
expressions, it is possible to take rh • h = 0, in which case the 
actual magnetization vanishes, but rh = §3 still defines the 
plane of the spiral. 



To determine rh, we must minimize Eq. (42 1. Let us first 
consider the special case c\ = 0, c 2 > 0. Then we may pre- 
sume rh lies in the plane spanned by z and h. Taking the angle 
of rh with the z axis as 9, and the angle of h with the z axis 
as 9h, the energy is 



E. 



FM 



= — c 2 cos — —— cos [9 — 
AE 

= -Acos[2(9 - 9 Q )} + const. 



0, 



where 



A 



-acos 



-4h 2 +4hcos26» ft , 

1 + 2h 2 cos 29 h 
VI + 4h 4 + 4h 2 cos 26 h 



(43) 
(44) 

(45) 
(46) 



with h = hj \J c-iAE. The angle #0 obviously gives the ori- 
entation of rh. Interestingly, it is an analytic function of h 
except at 9h = tt/2, i.e. when the magnetic field is per- 
pendicular to the (100) axis. As this value of 9h is ap- 
proached, #o(h) becomes sharper and approaches a step func- 
tion: 9 (h;9 h = n/2) = |0(h - 1/V2). It is also instructive 
to plot the magnitude of the magnetization, m(h). The mag- 
netization jumps at h = 1/ \/2 for 9h = tt/2, but is otherwise 
continuous (see Fig.[6]l. 

Before ending this section, we comment on the range of va- 
lidity of the results. First, though we have assumed through- 
out the above that c\ — 0, in fact it is possible to show that 
a spin flop (discontinuous jump in the magnetization) occurs 
throughout region I of the phase diagram in Fig. [3] in which 
\c\\ < c 2 . For brevity, we do not give the (algebraically in- 
volved) argument here. Second, we have assumed a particular 
ordering wavevector along the (110) direction. At zero field, 
this wavevector is chosen spontaneously from among the fam- 
ily of equivalent (1 10) planes (e.g. (101) etc.). In the presence 
of a field, the different wavevectors will become inequiva- 
lent, due to the magnetic anisotropy terms in /no- Hence, 
given enough time, annealing, or field cycling, the system may 
choose the lowest free energy wavevector amongst this set in 
the presence of the field. This is rather clearly the wavevector 
which is closest to the field axis. In this situation, the situation 
9h = tt/2 is avoided and the spin flop is avoided. In practice, 
wavevector reorientation is probably sufficiently slow at low 
temperature to allow observation of the spin flop transition. 
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FIG. 5: Ground state angle 8q of the magnetization versus field h for 
9 h = 0, 0.17T, 0.27T, 0.37T, 0.47T, 0.4997T (from the bottom curve to 
the top curve). 



m 



1.4 




- 0o=O 




1.2 




-- O.lff 

- 0.27T 




1.0 




-- 0.3tt 




O.S 




0.4tt 




0.6 




— 0.499^ 




0.4 


- 


0.2 


^> 






0.2 0.4 




0.6 



O.S 



1.2 



1.4 



FIG. 6: Magnitude of the magnetization m versus dimensionless 
field h for h = 0, O.Itt, 0.2tt, 0.3tt, 0.4tt, 0.499tt (from the top 
curve to the bottom curve). 



IV. COMMENSURABILITY EFFECTS 

Up to this point, our phenomenological theory leaves the 
phase of the spiral (i.e. the phase of d) free. In general, the 
different directions within the spiral plane are not equivalent, 
and when a fully account is taken of spin anisotropy and crys- 
tal symmetry, the phase of the spiral may take preferred val- 
ues. In this section, we discuss the effects of "pinning" of 
the phase, and how this leads to a "lock-in" transition for the 
spiral wavevector in some situations. 

We will assume the spiral form in Eq. Q, with some given 
q(J2/Ji) chosen to minimize the energy of the Heisenberg 
Hamiltonian. Using the arguments in the previous sections, 
we can fix the plane of the spiral, defined by the normal vector 
e 3 . Choosing two arbitrary unit vectors spanning the plane 
(§i x 62 = §3), we can then write 



d = m s (ei + ie2)e l 



(47) 



The terms considered up to now do not fix the phase 8. 

The freedom to choose 8 is related to translational invari- 
ance. In particular, under a translation r — * r + a, we have 



Here a can be any Bravais lattice vector. It is sufficient to 
consider the primitive lattice vectors a = (0, |, |) and per- 
mutations. We would like to construct terms in the effective 
continuum Hamiltonian or Landau free energy that are invari- 



ant under Eq. (48 I, but which depend upon 8 directly and not 
only through its gradients. Moreover, they must also be peri- 
odic in 8 (since a shift by 2tt leaves d unchanged). A general 
periodic functional of 8 can be written 



00 . 

Ve = -^2 d 3 r\ n cos(n8 + c/) n (r)), 



(49) 



where the X n are arbitrary coefficients. The <fr n are arbi- 
trary slowly varying functions of r, which should be chosen, 
if possible, to ensure invariance under Eq. (48 1. The func- 



tions should be slowly varying because large gradients of 8 
are heavily penalized by the Heisenberg Hamiltonian, which 
favors constant 8. If 8 varies slowly but <fi n varies rapidly, then 
this term will average rapidly to zero on integration, and can 
be neglected. 

A general choice of function which achieves the desired in- 
variance is 4> n (r) = nq • r(mod27r) + <fi n0 , with <fr n0 a con- 
stant. We need to determined which (if any) of these functions 
is slowly varying. This occurs if the change in </> n (r) on shift- 
ing by a primitive lattice vector is small. By continuity, this 
is achieved when q is close to a wavevector for which cf> n (r) 
is constant under such a shift. To achieve constancy, the n th 
term should have nq • a a multiple of 2tt for all three primi- 
tive vectors a. For this condition to hold for any n, we require 
that q • a be a rational multiple of 2ir. We call these special 
wavevectors satisfying this condition commensurate. 

Let us now specialize to a specific direction of wavevector 
of interest. We take q = (q, q, 0), corresponding to J2/J1 > 
0.7, which is the case appropriate for MnSc2S4. In this case, 
for the three primitive translations, we have qa = q/2, q/2, q. 
Thus the condition for the wavevector to be commensurate is 
nq/2 = 2wm, where n, m are integers. We assume the system 
is close to such a value, i.e. 



4 arccos 




4nm 



(50) 



q-a. 



(48) 



In general, the most important m, n will be those with the 
smallest n, since the terms A„ may be expected to decay with 
increasing n. For J2/J1 > 0.7 such that the (q,q,0) order 
is obtained, we find a number of commensurate wavevectors, 
shown in Fig. [2] The smallest q in this set is q = q^ s = 3tt/2, 
which is the wavevector observed in MnSc2S4 . The presence 
of these other commensurate wavevectors with smaller n sug- 
gests that other commensurate states might well be found by 
varying J2/J1 by physical or chemical pressure. 

Let us fix on the vincinity of one of these wavevectors. Be- 
cause the other terms in Vq rapidly oscillate, we need only 
keep the one involving q mn - 



(51) 



Vg = —A / d r cos(n8 + nSq ■ r + 4> n o), 



8 



where Sq = q q m ,„, with q m ,„ = (<?,„,„, q m ,n, 0), and we 
simplified A„ — > A. This term favors configurations in which 
V6> = (5<7 — 0„ o , which minimize the cosine. Establishment 
of a phase gradient, however, costs exchange energy. This can 
be seen because from Eqs. 1 47|2 1, a non- vanishing gradient 
V6> corresponds to a shift of wavevector. Indeed, the physical 
wavevector k for general 9 is 



k = q + V6 



(52) 



The exchange energy cost to distort the wavevector from q to 
k is, from Eq. ([Vj 



where the tensor stiffness is 




(54) 



Here k% = n + + k_, k 2 = k + — k_, and K3 are the stiffnesses 
along the principal axes. At zero temperature, they are given 
by: 



Ki = 



K-2 



Jl JJ 

2 16J 2 
(8J 2 - Ji)JiJ 3 



0[J 3 



(55) 



16Jf 



«3 



(128J| - 112JiJ| 



20 J X 2 J 2 - Jj 5 ) 



&J1J2 



For the most interesting case 5 « 37r/2, we have ki (2 + 
\/2) J1/8, k 2 m J3/2, and k 3 « (%/2 + 1) J 3 . 

We now proceed to analyze the effective Hamiltonian 
H e ff — H ex + Vg. Though we have given these expressions 
explicitly at T = 0, the general form in Eqs. ( 51|53 1 holds 
at any temperature below the Neel temperature, with H e ff 
replaced by F e tt, the effective free energy, and with renor- 
malized parameters A(T),k.;(T). Moreover, because in this 
temperature range the system exhibits magnetic long range 
order, the fluctuations of 9 are small and bounded, so that it is 
sufficient to consider saddle points of the free energy. 

It is convenient to shift variables to 9 = 9 + dq ■ r + <j> n Q jn. 
The free energy is 

F eff = [(fr^aJdJ + S-Vd 

-Xcosn9 +^Sq fl K fl „Sq u } 7 (56) 

with 8^ = K^Sq^. The last term is independent of 9 and can 
be neglected. One can readily see that 



which is along the (110) axis. We therefore define the coor- 
dinate x = (x + y)/ \/2 along the (110) direction, and rewrite 
the free energy accordingly, 

F eff = A J dx^idj)) 2 + 5dJ)- \cosn9}, (58) 

where A is the area of the sample transverse to the (110) axis, 
5 = \6\ = q — q m .n, and we have dropped the constant term. 
It is now evident that 5 enters only as a boundary term, which 
means that the free energy depends upon S only through the 
"winding" number N w = [6>(x = L) - 9(x = 0)]^ of the 
minimum energy saddle point (across the length L along the 
(110) direction). This allows one to proceed by finding the 
saddle point energy for fixed N w , and then minimizing over 
N w . 

It is useful to consider the cases N w = and N w = ±1. 
For N w = 0, the saddle point is uniform 9 = (up to a 
multiple of 2tt jn). For N w = ±1, one has a single soliton 
solution: 



9(x) = — arctan 
n 



±n 



(x-x ) 



(59) 



where Xo is arbitrary and specifies the location of the center 
of the soliton. Note that the soliton width w = — \f^f. The 

n V A 

energy of this solution, for 8 = 0, is EN m =i — -EW ra =o = 
8y/KiX/n. When the spacing between solitons is much larger 
than w, i.e. L/\N W \ 3> w, the energy of an N w soliton state 
is approximately just | | times this single soliton energy. 
Corrections to this non-interacting soliton approximation arise 
due to the overlaps of the exponential tails of the solitons. 
Defining the mean soliton density as n w = N w /L, we may 
then write the free energy density as 



2ttS 



-n w + c\n v 



(60) 



n n 

where c is a positive constant. From Eq. ( |60"|l, th e minimum 
n w can be easily found. For |<5| < |<5 C | = A^/kiX/tt, one has 
n w = 0, and the wavevector is commensurate. For \S\ > \5 C \, 
n w ^ 0, and the wavevector becomes incommensurate. Due 
to fluctuations, one expects both A and Ki to decrease with 
temperature. Hence the width of the commensurate state (cx 
\S C \) will decrease with increasing temperature. A schematic 
phase diagram is shown in Fig. [7] 

We see that, when J2/J1 is close (but not too close!) to 
a value for which the Heisenberg model alone has a com- 
mensurate spiral solution, there is a "lock-in" transition on 
decreasing temperature from an incommensurate to a com- 
mensurate spiral. Within the commensurate ("C" in Fig. [7]) 
phase, the wavevector is constant and equal to q myn . This is 
consistent with observations on MnSc2S4 . Commensurate- 
Incommensurate transitions of this type are well studied, and 
the reader interested in details of the associated critical behav- 
ior may find it in various standard texts, for instance Ref.[7j 



k = q„ 



(57) 



V. QUANTUM FLUCTUATIONS 



The minimum free energy saddle points of F e ff are trans- 
lationally invariant along the directions perpendicular to d, 



In this section, we develop a spin wave theory for the di- 
amond antiferromagnet, and obtain the leading quantum cor- 



9 



B. Spin wave Hamiltonian 




Inserting Eq. (j64b into the Heisenberg Hamiltonian, we ob- 
tain terms of ©(S 2 ), <D(S 3 / 2 ), and O(S), dropping higher or- 
der corrections: 



H. 



o(s 2 ) 



— ~ S' J' iiZj, ' Z 



(66) 



o(s 3 / 2 ) 



2V2 



H (S) = ■^SJ lJ [{n l + n, J )i l -i, J +x l X]^i-^j+^ij\- 

The 0{S Z / 2 ) term vanishes because the local coordinate 
vectors are eigenstates of the exchange matrix, e.g. 



FIG. 7: Schematic phase diagram showing commensurate (C) and 
incommensurate (IC) magnetic phases, and the paramagnetic (PM) 
phase. The figure is drawn as appropriate for a first order magnetic 
transition line, in which case the width of the commensurate phase 
remains non-zero on approaching the Neel temperature. 



rections to the spin correlations. 



A. Holstein-Primakoff bosons 



(67) 



and Xi • ii = 0. Here J m is the minimum eigenvalue of 
the exchange matrix. The vanishing of the 0(S 3 ^ 2 ) term is 
of course true because we expand about the classical ground 
state. 



Using Eq. ( 67 1, one can further simplify the spin wave 
Hamiltonian. We obtain 



H, 



O(S) 



(68) 



■ E ^r(x? + a + E ■ + ^ - 



We proceed in the standard way by defining Holstein- 
Primakoff bosons in a spin coordinate frame rotated to follow 
the classical ordered state. The local orthonormal axes will be 
defined by 



X; 



Sf = Re [de* q ' r * 
-Im[de**'*], 



(61) 



d x d* 

2 



e 3 - 



Note that the y axis is site independent, as it just corre- 
sponds to the normal vector to the spiral plane. The linearized 
Holstein-Primakoff transformation is 



{S- ni )zi + V2S(al 



(5q + iy) 



(x; - iy) 



), (62) 



which neglects corrections cubic in the canonical a,{ , a] bo- 
son operators. Here m — ajc^ as usual. It is convenient to 
pass from canonical bosons to "coordinate" and "momentum" 
operators, 



1 / t^ 



£ i = i—={a\-a,i). (63) 
v 2 



The spin operator becomes 

Si = (S - n^Zi + (xA + &y) 

and 

& 1 



n. 



2 



(64) 



(65) 



neglecting constant terms which do not affect the correlations. 



C. Action 

Spin fluctuations are conveniently calculated using the path 
integral approach. The imaginary time action corresponding 



to Eq. ( 68 1 has the usual Berry phase terms describing the 
canonical commutation relations of Xi>£i> 



S 



J j#o(S) +E% 9 t& j 



(69) 



Static correlations of Xi ar, d £?' vanish, so we may consider 
the two separately. It is then convenient to integrate out one 
of these fields to obtain an effective action for the other. This 
gives 



S-y — 



S f = 



^SlJ {SJiMi + ^[K-%d T ^}, (71) 



where 



J^j J ij J m Si j , 



(72) 
(73) 
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To diagonalize this, we move to momentum space. Due to 
the sublattice structure, we define two components for each 
field, XAk, Xsk and £_ Ak , £ Bk , such that 



L 



Xi = I Xs(i)ke lk ' r \ 

'k 



(74) 
(75) 



with s(i) — A,B specifies the diamond sublattice of the site 
i. The k integral is defined as J k = v uc J £J? 3 , where the 
integration domain is the first Brillouin zone, and v uc = 1/4 
is the volume of the real space unit cell. It is convenient to 
define 



Xk 



XAk 
XBk 



£sk 



(76) 
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FIG. 8: Reduced magnetic moment M s as a function of J3/I Ji|. 
The solid line is for fixed J%, while the dashed line is for J 2 , Jz satis- 
fying Eq.[TT] so that the wave vector remains equal to 37r/2(l, 1, 0). 



The action becomes 

x 2 
1 



X 

-1/ 



(77) 
(78) 



Here the frequency integral is J = J ~ as usual. The matrix 
Green's functions are straightforwardly found, but somewhat 
cumbersome. The reader interested in the details is referred 
to Appendix [B] With all these definitions, one can formally 
evaluate the equal time correlation functions: 



(XiXj) 



[Gx(k^)L 0>w e lk -^\ (79) 



[Gf(k,u;)] 



s(j)s(i) 



ik- (r^ — r j ) 



Here the subscripts give the matrix elements of the matrix 
Green's functions. 



D. Local moment 

Focusing on the case of M11SC2S4 , with q = 
(3tt/2, 37t/2, 0), we have calculated the reduction of the sub- 
lattice magnetization by numerically evaluating the momen- 



tum integrals in Eqs. (79i (the frequency integration can be 
done analytically). See Appendix [C] for more details of the 
calculation. The result for the on-site expectation value is 



(X? 



0.67, 



(3 



1.19, 



(80) 



for J 3 = 0.1K ps Ji/100. From this, one obtains {n % ) pa 0.43 



from Eq. (65 I, which is approximately a 20% reduction from 
the classical local moment. As J 3 is increased, the moment 
increases closer to the classical value, as shown in Fig. [8] 



VI. 



MICROSCOPIC ORIGIN OF MAGNETIC 
ANISOTROPY 



In Sections II B III we studied the effects of explicit spin 
rotation symmetry breaking on phenomenological grounds, 
using only the space group symmetry of spinel structure. In 
this section, we address its microscopic origins. There are 
in general two mechanisms of spin rotation symmetry viola- 
tion in solids: (1) dipole interactions between electron spins, 
and (2) spin-orbit coupling. We consider both in turn, and 
find these lead to somewhat different regimes of the phe- 
nomenological model discussed previous. Interestingly, only 
the spin-orbit coupling mechanism can explain the observa- 
tions in MnSc 2 S4 . 



A. Dipolar interactions 

The dipole-dipole interaction can be written as 
/in m, • Hi, 



H, 



3mj • Tijiaj ■ Yij 



4ir ^ 

1,3 



(81) 



where = g/isSSi is the dipole moment of the spin i (we 
included an explicit factor of S to follow our convention of 
unit vector spins). Using g ps 2 as expected for an S = 5/2 
Mn 2+ spin with a half-filled d shell, we obtain a dipolar en- 
ergy of interaction between two nearest-neighbor spins of ap- 
proximately 0.5K. We note that this is not negligible (espe- 
cially when added over many spins within a correlation vol- 
ume) but it is certainly weak compared to the basic energy 
scale of exchange interactions as estimated from the Curie- 
Weiss temperature Qcw ~ — 23K. Therefore we expect 
we can treat the dipolar interaction as a weak (but symme- 
try breaking) perturbation on the ordered ground states of the 
Heisenberg model. 

To this end, we first consider the dipolar interaction classi- 



cally by simply inserting the general spiral form of Eq. ( 29 1 
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in Eq. ( 8 1 1 and evaluating the sum. Because we are only in- 
terested in the dependence of the energy upon the spin orien- 



tation of the spiral, we may drop the first term in Eq. (81 
which is fully SU(2) invariant. Because the spiral itself is at 
a non-zero wavevector, there are no convergence difficulties 
with the long-range dipolar sum. Choosing the wavevector 
q = (q, q, 0) as in experiment, one indeed finds the form in 



Eq. ( 22 1 is obtained provided the sum is truncated in a manner 
preserving cubic symmetry. We plot the values of c\ and C2 
in Eq. p7| ) in the physical range of q for 0.7 < J%j J\ < 1 in 
Fig. [9] Throughout this range we find c\ > and more than 3 
times as large as c%. This favors alignment of spins within the 
plane normal to 63 = (110). Unfortunately, this is not what is 
found experimentally. 



dipole-dipole interaction as a perturbation, it is sufficient to 
consider the expectation value (Hp) in each of the spin wave 



ground states specified by d. To do so, we insert Eq. ( 64 1 
into He,, and expand to quadratic order in Xi>£i> then take the 
expectation value of the result. The necessary correlators of 



Xi, £i are calculated by numerical integration of Eq. (79 1. The 
values obtained are given in Table[I] Because the basis vectors 
Xj, y, Zi are expressed in terms of d in Eq. (61 1, the result is 
again an energy function of the form of Eq. ( 27 1, which con- 



tains both the classical expressions for ci , C2 and their leading 
quantum corrections. We find that the quantum corrections 
push the system even further from the 63 = (110) state, and 
in any case the magnitude of the corrections are very small 
compared to the classical values. 



Cl,C 2 

-75- — 
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(XiXj) 




0.67 


1.19 


0.22 


-0.1 


0.18 


-0.3 


-0.25 


0.23 


-0.25 


0.23 





±1(1,1,0) 
§(±1,0, ±1) 

1(1,1.-1) 
-£(1,1,1) 



TABLE I: Numerically calculated values of correlations of f»,Xi 
fields from spin wave theory, for q = (3tt/2, 3n/2, 0), J2/J1 = 
! cos 2 (7r/8), J3/J1 = 0.01. Values not specified have negligible 
correlations. 



0.70 



0.72 



0.74 



0.76 



FIG. 9: Calculated anisotropy parameters (in arbitrary units) ci (up- 
per curve) and C2 for q = (q, q, 0) as a function of qj (27r). 

Several possible complications should be considered before 
abandoning dipolar interactions as a mechanism of magnetic 
anisotropy. First, in applying Eq. (81 1 with m, = gfigSSi, 
we have treated the electron spins as point dipoles. In fact, 
the electronic wavefunctions may be somewhat extended. 
Through such "covalency", there may be some spin density 
not only in the atomic d orbital of the Mn 2+ ion, but also on 
the neighboring chalcogenide p orbitals. This can be approx- 
imately accounted for by modifying the dipole moment dis- 
tribution associated with a spin accordingly, to be distributed 
amongst with a fractional moment 1 — / on the central Mn 2+ 
ion and a fraction //4 on each of the neighboring four S 2 ~ 
ions. We have carried out such a modified dipolar sum, and 
found that it does not substantially alter the results of the point 
dipole model for a reasonable range of parameters /. 

Another more interesting possibility is that fluctuations 
may alter the dipolar energetics. This is not an unreason- 
able possibility to consider since, although the classical or- 
der parameter description is expected to qualitatively (and in- 
deed rather quantitatively) capture the long-range order of the 
spins, the dipolar energy actually receives large contributions 
from very nearby spins. The latter could exhibit quite differ- 
ent correlations from well separated spins which control the 
order parameter. 

To consider this effect, we have calculated the leading cor- 
rections in 1/5 to the dipolar energy using the spin wave for- 
malism described in the previous chapter. Since we treat the 



Having thus exhausted the possible complications associ- 
ated with the dipolar interactions, we conclude that the ob- 
served ordered state in MnSc2S4 is inconsistent with a dipolar 
origin of the magnetic anisotropy. We therefore turn to spin- 
orbit effects in the following subsection. 



B. Exchange anisotropy due to spin-orbit coupling 

As we saw in the previous subsection, dipolar interactions 
do not appear to be viable explanation of the orientation of 
the spin spiral observed in MnSc2S4 . We now consider the 
second microscopic origin of magnetic anisotropy, which is 
spin-orbit coupling. From the point of view of symmetry, the 
spinel lattice allows both single-ion (cubic) anisotropy of the 
Mn 2+ spins and exchange anisotropy. The former is how- 
ever expected to be extremely small for Mn 2+ , which has an 
extremely stable and isotropic 3d 5 configuration (one may ex- 
pect a coupling constant of a few millikelvin). However, ex- 
change anisotropy is non-negligible in many Mn magnets. A 
microscopic calculation is beyond the scope of this paper, but 
we can make a few statements on general grounds. Because 
of the closed shell configuration, these effects are also ex- 
pected to be much smaller than the typical exchange interac- 
tions (i.e. perturbative in spin-orbit coupling). However, they 
may still be as large as or larger than the dipolar effects. In 
MnSc2S4 , one may attempt to get some feeling for their mag- 
nitude by comparing the measured effective moment seen in 
the Curie law /i e f t = 5.8/ib to the theoretical spin-only value 

^s=b/2 = 2y^f| ~ 5.92. Given uncertainties in the mea- 
surement, we expect no more than a 5 — 10% deviation from 
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the latter (and very possibly a much better agreement masked 
by experimental complications). For Mn 2+ , one expects that 
contributions to the g-factor (which renormalize the effective 
moment) are second order in the spin-orbit coupling. Ex- 
change anisotropy occurs at both first order and second order. 
At first order, one obtains the antisymmetric Dzyaloshinskii- 
Moriya (DM) interaction, and at second order symmetric ex- 
change anisotropy. Thus we would expect that the DM inter- 
actions be of order y/\^ e ff — M5/2 1 / ^5/2 Jij an d symmetric 
exchange anisotropy be of order (|/i e // — M5/2 \/ 1^5/2) Jij • 

With this in mind, we consider the allowed form of the ex- 
change anisotropy as constrained by the space group symme- 
try of the spinel structure. We first consider nearest-neighbor 
bonds. Without loss of generality, take a bond oriented along 
the (111) axis. DM interaction is forbidden on this bond be- 
cause exactly between the two sites is an inversion center (Q§ 



in Eq. ( 17 1). Thus we need only consider exchange anisotropy. 
This in turn is strongly constrained by the C3 rotation sym- 
metry about the (111) axis {Qj, in Eq. ( 17 1). This allows only 
two separate exchange couplings, for components parallel and 
perpendicular to the bond, respectively. We can write the as- 
sociated exchange Hamiltonian as 



rrn.n 

ani 



J±ja.n x Si 



(hi) 



(82) 



There is a single parameter, J± — J», which parametrizes the 
nearest neighbor exchange anisotropy. 

Next, we consider the anisotropic exchange for next- 
nearest-neighbors. Here the symmetry is considerably less 
constraining, since two second neighbors (fee neighbors) are 
not connected by a C3 axis, and there is no inversion center 
between them. We have however determined the most general 
exchange Hamiltonian between two such sites invariant un- 
der all operations in Eq. ( |T7) , which is a straightforward but 
tedious calculation. There is unfortunately no simple expres- 
sion for this Hamiltonian which describes all 6 second neigh- 
bor bonds simultaneously. Instead we write the form for a 
particular bond, connecting two sites i, j on the "A" sublat- 
tice, separated by the (arbitrarily chosen) Bravais lattice vec- 
tor Tij = (0,-i, |): 

rrn.n. n j qx qx , t-\( qx qV q%) qx , qx qz qz qx\ 

+J b (SfS? + S*S]) + J c (SfS] + SIS*). (83) 

The full set of Hij for all other pairs of second neighbor 
sites can be obtained by actions of symmetry operations on 
Eq. ( [83) , which thus defines the full next nearest neighbor 
Hamiltonian H a \™ n -. Note that there are 3 symmetric ex- 
change constants, one linear combination of which represents 
the isotropic Heisenberg term, and the other two ( J , J a — J c ) 
represent symmetric exchange anisotropy. Because of the 
absence of an inversion center between two fee sites in the 
spinel, there is an allowed DM term D. However, the pres- 
ence of the inversion center implies that the D term takes the 
opposite sign for spins on the "B" sublattice. 

We can now consider the full exchange anisotropy Hamil- 
tonian, Hani = Hani + H ani' n ~ #Heis., as a perturbation 



to the Heisenberg form, and evaluate the energy splittings in- 
duced for a given spiral state specified by q and d, by sim- 
ply inserting Eq. Q in -ff an ;. As required by symmetry, for 
q = (qqO) it again has the form of Eq. (22 1. Reading off the 
coupling constants, we find 



ci = (J a -J c )(l + V2) + {1-^=)(J ± -J l{ ), 



C-2 



Jb- 



(84) 



Note that the DM term D does not enter these macroscopic 
anisotropy parameters, which is a consequence of its stag- 
gered nature on the two diamond sublattices. 



Unlike for the dipolar interactions, we see that Eqs. (84i 
allow essentially arbitrary values of c\ and C2. This means 
that, in the absence of a microscopic calculation, the exchange 
anisotropy mechanism is not inconsistent with the observed 
ordering in MnSc2S4 , which as we saw could be described 
phenomenologically by a range of choices of c\ , C2. Given the 
incompatibility of our dipolar results, however, we tentatively 
conclude that spin-orbit induced exchange anisotropy is likely 
at the origin of spin state selection in MnSc 2 S4 . 



VII. DISCUSSION 
A. Summary 

In this paper, we have extended the theory of Ref. 1 to de- 
scribe the effects of magnetic anisotropy and quantum fluctu- 
ations in frustrated antiferromagnetic A-site spinels. The the- 
ory predicts the possible planes in which spins reside in the 
spiral magnetic ground states in zero field, and describes their 
evolution with field. In some orientations a spin flop transition 
was found. We described commensurate-incommensurate 
transitions which occur below the Neel temperature when the 
spiral wavevector locks to one of a set of specific commen- 
surate values. These effects are all in accord with observa- 
tions on the best studied such material, MnSc2S4 . We ad- 
dressed the reduced static moment seen in MnSc2S4by spin 
wave calculations, and found that a relatively large reduction 
can indeed be achieved by quantum fluctuations due to the 
frustration-induced degeneracy, despite the large 5 = 5/2 
spin of Mn 2+ , if one assumes the third neighbor exchange 
J3 < 0.1. Finally, we derived microscopic expressions for the 
most important phenomenological magnetic anisotropy pa- 
rameters, taking into account both dipole-dipole interactions 
and spin-orbit effects. In MnSc2S4 , we concluded the latter 
are most likely responsible for the observed magnetic orienta- 
tion. 



B. Experiments 

Let us turn now to a further discussion of experiments. First 
we discuss existing results, and then consider future experi- 
ments. 
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1. local moment 

As mentioned above, from the weight in the magnetic 
Bragg peaks seen in Ref. 6 in M11SC2S4 , it was estimated 
that the local ordered moment M s w 0.8M C ;, where M c i is 
the expected classical static moment for an S — 5/2 spin. 
In Sec. VD we showed that the 17% reduction could per- 
haps be due to quantum fluctuations, if J3 is sufficiently small. 
However, there are a number of reasons to be cautious about 
this conclusion. First, at a technical level, it is not clear to 
us how large the experimental errors should be considered on 
this measurement, which was done in a powder sample. Sec- 
ond, the data was taken at T — 1.5K, more than half the 
ordering temperature T c = 2.3K, so thermal fluctuations may 
contribute to some reduction of the moment. 

Finally, there are a number of different effects that have not 
been addressed theoretically, which may contribute to the mo- 
ment reduction. First, we have neglected disorder, which is 
known to be present in the form of inversion - interchange of 
A and B site atoms of the spinel. Such disorder can damage 
the spin spiral, reducing the ordered moment even if the lo- 
cal static moments remains large. The nature of the defects 
created and their impact on the ordered moment measured by 
neutrons will be discussed in a separate future workP A sec- 
ond effect that could contribute is a spin-orbit renormalization 
of the g-factor. Usually this is small in Mn 2+ magnets, but 
perhaps this is something worth considering further. 



2. microscopies ofanisotropy 



VI though dipolar inter- 



As discussed above and in Sec. 
actions between Mn 2+ spins might seem a likely candidate 
for the origin of the magnetic anisotropy in MnSc2S4 , they 
appear to be inconsistent with the observed nature of this 
anisotropy. While we can reconcile the existing experiments 
with a picture of spin-orbit induced anisotropy (with some as- 
sumptions), it is still surprising to us that such effects would 
be competitive with dipolar interactions. We believe the con- 
flict of the latter with the ordered state seen in MnSc2S4 is a 
significant one, and found in Sec. | VI A| that neither covalency 
nor quantum fluctuations were likely to effect a reconciliation. 

One possibility we have not considered is the effect of dis- 
order and granularity. Given the long-range nature of the dipo- 
lar interaction, it is possible that defects created by disorder in 
an ideal spiral can facilitate large changes in the dipolar en- 
ergy. This is an interesting issue to be explored in the future. 
We emphasize that, although such a mechanism of anisotropy 
might be possible, the phenomenological portion of our theory 
is entirely independent of these details and is quite generally 
valid irrespective of the microscopic physics of anisotropy. 



3. magnetization experiments 

We now turn to future experiments. Of particular interest 
would be the development of single crystals. This has already 
been emphasized in Ref. Q] where predictions were made for 



unusual "spiral surface" structure in the angle-resolved neu- 
tron structure factor. Based on the results of this paper, we 
suggest that single crystals are also interesting for the study 
of magnetization effects. An obvious suggestion is to look for 



signs of the spin flop transition discussed in Sec. IIIB An- 
other interesting measurement would be torque magnetome- 
try. As shown in Fig. [5] the angle of the magnetization can 
be strongly misaligned with the applied field, which should 
leads to a large torque. This is a very sensitive technique that 
perhaps does not require as large crystals as neutron scattering 
does. 



C. Ferroelectricity 

Our results enable us to discuss magnetically-induced fer- 
roelectricity in the A-site spinels. This may be expected 
since many recent studies, both theoretical and experimen- 
tal, have emphasized the relation between spiral spin states 
and ferroelectricity. The basis for such a relationship goes 
back much earlier to symmetry considerations of Landau^ and 
Dzyaloshinskii.L^ Several recent studies have pointed out that 
very general arguments suggest a simple relationship between 
the electric polarizat ion P a nd the basic parameters 63 and q 
describing the spiralPDE] 



P oc e 3 x q. 



(85) 



Here, as in the text, §3 is the axis which is perpendicular to 
the plane of the spins and q is the wavevector. 



The argument leading to Eq. ( 85 1 is rather simplified, and 
actually assumes a sort of "spherical symmetry". In reality, 
in the reduced crystal symmetry environment of the solid, 
the actual relation may be somewhat different. Still, for the 
A-site spinels, a complete symmetry analysis leads to rather 
similar results. In particular, time reversal symmetry allows 
a quadratic term in the d order parameter (which is time re- 
versal odd) to couple linearly to P. One therefore expects the 
polarization to take the form 



P a = c Q ^ 7 (q)d^d 7 . 



(86) 



As argued earlier, all such bilinears in d can be rewritten in 
terms of §3. The coefficients c Q ^ 7 are constrained by crystal 
symmetry. Specifically, we require that the left and right hand 
sides of Eq. < [86] > transform identically under the little group 
which leaves q invariant. 

For q = (q, q, q), applying Eqs. ( 19 1, we find the form 



P111 = ci 




(87) 



The simplified Eq. ( |85) corresponds to ci = 0, C3 = — c 2 . 
However, in general, symmetry allows any values of ci, C2, C3. 
For q = (q, q, 0), using Eqs. (20i, we find instead 



P110 = c i 



C2 



(88) 
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Eq. (85 1 is the special case c-i = — c\. 

Given these results, we can make some limited predic- 
tions on the ferroelectric polarization in the A-site spinels. In 
MnSc2S4 , where the ordering wavevector and spiral plane is 



known, we can directly apply Eq. ( 88 i without much ambigu- 
ity. We have §3 = z, which means that there is a spontaneous 
polarization with P along the 110 direction. It would be inter- 
esting to search for this experimentally in single crystals, or 
for dielectric anomalies related to this in powders. Moreover, 



the phenomenological theory in Sec. Ill B in conjunction with 



Eq. (88 1, describes how this polarization may be rotated by an 



applied field. Again, detailed single crystal studies would be 
enlightening. 

For spinels in the regime where q = (q, q, q), the theory 
is somewhat less predictive. This is because not only is there 
ambiguity in the spiral plane giving 63 (due to the unknown 
constant c in Eq. (26i), but also there are more unknowns 



in the relation between the polarization and the spiral plane, 
Eq. (87 1. A microscopic theory for Eq. (87i, which deter- 
mines the Ci, is therefore desirable. We imagine one might be 
constructed based on the inverse Dzyaloshinskii-Moriya inter- 
action mechanism 12 , since we have seen that there is a single 
DM interaction allowed in the A-site spinels - see Eq. ( S3 1. 



The polarization can be very sensitive to details of the mi- 
croscopies. For instance, for c > in Eq. (26 1, we have 
63 = (1, 1, and according to Eq. (85 1, the polariza- 

tion vanishes. However, in general this is an artifact of the 
simplifications in Eq. ( 85 1, and according to Eq. ( 87 1, P ^ 0. 



However the orientation of the polarization is precisely con- 



trolled by deviations from the naive Eq. ( 85 1. 
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value of k z into Eq. ( |10) l, we can obtain the energy on the sur- 
face explicitly. One can then scan linearly along lines defined 
by k = (qcos9, qsin9, k) on the surface and determine the 
lowest energy for each 9. In every case, the lowest energy as 
a function of 9 is achieved for 9 = tt/4 (an example is shown 
in Fig. 10 1, which implies a wavevector of the form (q, q, k) 
on the surface. 




FIG. 10: Minimum value of S(k) for k of the form k = 
(qcos8, qsin9, k) as a function of 9, for J2/J1 = 0.4. For all 
values of J2/J1, the minimum value is achieved at 9 — 7r/4. 

Having determined that the ground state wavevector is al- 
ways of the (q, q, k) form, we need only search this ray for the 
ground state. This can be done analytically. One obtains 

= Z(q,q,k z (q)) 
' |96(8A 2 - 3) A 2 + (256A 4 - 13) cos q 



64 

+2(16A 2 4- 5) cos 2q - 3 cos 3q + 6 



sec 2 -. (A2) 
2 



This should be evaluated only when such a wavevector exists 
on the surface. This condition is 



sin 2 \ < A 
4 



sin 2 q - > 1 
4 



A. 



(A3) 



APPENDIX A: SPLITTING OF SPIRAL SURFACE 
DEGENERACY 

In this appendix, we give some details on how the ground 
state spirals are determined in the presence of third neigh- 
bor antiferromagnetic exchange J 3 . First, we performed a 



numerical study of the minima of Eq. ( 10 1, considering only 
wavevectors fixed on the spiral surface, i.e. satisfying A(k) = 
A = 1/8 J2. These can be conveniently studied by solving this 
condition to give k z in terms of k x ,k y : 



±4 arccos 



A 2 



sin 2 'f sin 2 'f 



v cos 2 ^ cos 2 % 



sin 2 f sin 2 k -f 



1/2 



(Al) 



Here the solution (and the surface) exists only when the ar- 
gument of the square root is between and 1. Inserting this 



Now it is simple to study the ground states. One can check 
that the (q, q, q) state, for which q = arccos[(8A 2 — 5)/3], is 
always a local minimum of Eq. (jA2j). It however only exists 
however when this value is well-defined, which requires A > 
1 /2. This corresponds to 1/4 < J2 < 1 /2. Indeed, in this 
range it is straightforward to show that this is the global energy 
minimum. 

For J2 sufficiently large, one can readily see that the min- 



imum of Eq. ( A2 1 is instead achieved at the boundary of its 
domain of validity, i.e. when the inequalities in Eq. ( |A3[ ) 
are satisfied as equalities. This corresponds to k z (q) = 0, 
i.e. a (q, q, 0) state. This eventually ceases to be a mini- 
mum for small enough J2- A choice of such wavevector is 
q = qo = 4 arcsin VA. For this to be a minimum, we need 
£'(<7o) < 0. By differentiating Eq. (A2i and evaluating, we 
find 



£'( 90 ) = 16 [A(l - A)] 



3/2 



1 - 6A + 4A 2 
2A- 1 



(A4) 
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It is straightforward to show that this is negative provided A < 
(3 - VE)/4 or J 2 > l/[2(3 - Vb)}, which determined the 
domain of the (q, q, 0) state. In between this and the (q, q, q) 
state we necessarily have the (q, q, q*) state. 



APPENDIX B: SPIN WAVE GREEN'S FUNCTIONS 



1 



1 



lu 2 1 



UJ 



+ PlUJ 2 +p 2 2y/pi + 2^/p2 



In this appendix, we give some details of the spin wave 
Green's functions. The Green's functions defined in Eq. ( |77] > 
can be written as 



The frequency integrated Green's functions, G x r^ (k), can 
be written as 



G x (k,w) = S , B(k) + w 2 (S'A(k))- 1 
G* e (k,w) = [s , t4(k) + w 2 (S' < B(k))- 1 



, (Bl) 



(B2) 



G 



x(€) 



(k) = 



(CI) 



with 



Here we've defined a number of matrices occurring as Fourier 
transforms of exchange matrices: 



G"(k) = G 22 (k) 



A(k) = W q , 7 (k), S(k) = W ,o(k). (B3) 

<-> / Wi 1 fk) Wj, 2 (k) \ 

W*(k) ee ( ^'rfj >L ) • (B4) 



C(k) 

21/ 



t^Jk) WW(k) 



The elements of are conveniently given in terms of the 
nearest-neighbor vectors n Q of the A sites of the diamond lat- 
tice, 



G' 2 (k) = [G 21 (k)] 



1 



C(k) 



(B 12 (k)-£>(k)A 12 (k)) 



(C2) 



(C3) 



(C4) 



n = |(1,1,1), ni = 1(1,-1,-1), (B5) Eq.[B3| Then, C(k),D(k) are, 
n 2 = i(-l,l,-l), n 3 =i(-l,-l l l). (B6) 



Here A "(k) is the a matrix element of A (k) defined in 



Then 



W k V 7 (k) = W# i7 (k) (B7) 
= -J™ + J 2 X] e ik ( n »- n ^ cosk' • (n Q - n 6 ), 



G(k) = 2^jTr[A(k)B(k)} + 2\A(k)B(k)\ 
D(k) 



S(k)| 



\ l^(k)l 



W^k) = [w^ i7 (k 



Ji^e*-"- cos(k'-n a + 7 ) 



(B8) 



+ ^ J 3 E e ik ' (n " +n ^ nc) cos[k'-(n Q + n fc -n c )+ 7 ]. 

Here the sums range over distinct values of a, 6, c taken from 

0,1,2,3. 



APPENDIX C: FREQUENCY INTEGRALS OF GREEN'S 
FUNCTIONS AND MOMENTUM INTEGRATION IN 
BRILLOUIN ZONE 

In this appendix, we give some details of the frequency inte- 
grals of Green's functions and the transformation to unit vari- 
ables Xi in momentum space. The frequency integrals of the 



It is natural from Eq. B 1 that G j (k) can be expressed 
G x (k) with the changes A (k) <-> B (k) 

The numerical integration of the momentum in the 1st bril- 
louin zone, can be eqsily evaluated using the transformation 
to the unit variables in momentum space. 

k = bixi + b 2 x 2 + h 3 X3 
_ 2-Ka.j x a fc 

* i \ 

ai ■ (SLj x an,) 

where is the fee primitive vectors, permutations of 
1/2(0,1,1). Hence we can transform the momentum k to unit 
variables x$, then the momentum integration in the 1st bril- 
louin zone can be written as 



correlation functions defined in Eq. ( 79 1 can be calculated an 



BZ 



d 3 k 

(2tt) 3 



n 



i=l 



dxi 



(C5) 



alytically using the following relations. 
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